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SUBMILLIMETER  WAVE  PRODUCTION  BY  UPSHIFTED  REFLECTION 
FROM  A MOVING  IONIZATION  FRONT 

Introduction 

An  ionization  wave  can  be  propagated  through  a gas  in  many  ways, 
e.g.  a strong  electromagnetic  pulse,  an  ionizing  shock,  a propagating 
electron  beam,  a laser  or  electron  beam  sweeping  across  the  gas,  or  a 
programmed  sequence  of  laser  pulses.1  The  velocity  U of  the  ionization 
front  can  be  relativistic,  or  even  (in  the  latter  cases)  exceed  c. 
Similarly,  a finite  time  after  such  an  ionizing  pulse  has  passed,  a 
recombination  wave  will  move  across  the  plasma.  In  this  paper,  we 
consider  the  reflection  of  an  electromagnetic  wave,  incident  from  the 
neutral  gas  side,  by  such  a moving  ionization  or  recombination  front. 

It  would  appear  that  an  electron  density  which  will  render  the  plasma 
overdense  to  microwaves  can  be  attained  with  relatively  modest  means. 

We  suggest,  therefore,  that  the  Doppler  shift  which  occurs  upon  re- 
flection can  be  used  to  upshift  microwaves  to  the  submillimeter  wave 
range.2 

For  specificity,  we  assume  that  the  ionization  or  recombination 
front  is  planar,  and  that  the  EM  wave  is  normally  incident.  In  the 
present  calculations,  the  plasma  is  treated  as  cold  and  collisionless. 
The  calculations  are  most  easily  performed  in  the  frame  of  reference  in 
which  the  ionization  or  recombination  front  is  stationary  [i.e.  the 
electron  density  n(z)  is  time- independent]  , and  the  gas/plasma  is 
streaming  at  velocity  U,  as  illustrated  in  Fig.  1. 

The  reflection  process  under  consideration  has  some  similarities 
to  and  some  differences  from  the  process  of  reflection  by  a moving 
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mirror.  A relevant  example  of  the  latter  is  upscattering  of  microwaves 


by  reflection  from  a relativistic  electron  beam.  In  both  cases,  re- 
flection is  from  the  electron  density  profile,  rather  than  from 
individual  particles,  and  thus  the  frequency  cur  of  the  reflected  wave 
is  Doppler  shifted  relative  to  the  incident 


tur/o)±  - (1  + 3)/(l  - 3), 


where  3 = U/c.  In  both  cases,  the  duration  of  the  reflected  pulse 
is  similarly  related  to  that  of  the  incident  pulse. 


VTi = ^ - 0)/(1  + *)• 


However  the  plasma  itself  is  moving  in  the  case  of  the  moving  mirror 
(e.g.  relativistic  electron  beam)  but  is  stationary  in  the  case  con- 
sidered here.  This  introduces  a difference  in  the  wave  dispersion 
properties  in  the  plasma  which,  in  the  case  of  an  oncoming  ionization 
front  (U  > 0 in  Fig.  1),  reduces  the  reflection  coefficient  as  compared 
to  what  it  would  be  for  an  oncoming  mirror.  Indeed  it  is  obvious  that 
this  must  be  true,  since  the  pulse  reflected  from  a moving  mirror  has 
more  energy  than  the  incoming  pulse,  the  extra  energy  being  supplied 
by  the  kinetic  energy  of  the  mirror;  in  the  case  of  reflection  from  a 
moving  ionization  front  in  a stationary  gas,  this  energy  source  is  not 
available.  However,  the  reflection  coefficient  from  a receding  recom- 
bination front  is  shown  to  be  exactly  the  same  as  that  from  a receding 


mirror. 
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Although  the  reflection  coefficient  is  smaller  for  an  oncoming 
ionization  front  than  for  an  oncoming  overdense  electron  beam,  the 
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efficiency,  i.e.  the  reflected  electromagnetic  wave  energy  divided  by 
the  total  energy  input  (incoming  wave,  plus  beam  energy  or  energy 
invested  in  ionization)  of  the  two  schemes  may,  in  practical  cases, 
be  comparable.  Furthermore,  the  ionization  front  scheme  has  advantages 
with  respect  to  size  and  cost. 

Analysis 

We  perform  the  analysis  in  the  frame  of  reference  in  which  the 
ionization  front  is  stationary  (Fig.  l).  We  make  no  assumptions  about 
the  shape  or  width  of  the  ionization  front;  it  may  be  sharp  or  broad, 
the  only  requirements  being  that  n(z)  -*  0 for  z - ®,  and  n(z)  const 
for  z -*  + ®.  We  let  ,j(z,t)  * - n(z)ev(z,t)  be  the  oscillating,  trans- 
verse plasma  current;  y(z,t)  is  thus  the  oscillating  component  of  the 
electron  fluid  (i.e.  local  mean)  velocity,  which  is  determined  by  the 
momentum  conservation  equation, 

6v  6v  dn  . 

n + nU  -g—  + uv  + - UxB ) + R,.  (3 ) 


v For  the  case  of  a recombination  front  (U  <0),  R = Uv(dn/dz)  repre- 

sents  the  loss  of  electron  oscillating  momentum  upon  recombination; 
r,  for  the  case  of  an  ionization  front  (U  > 0),  R = 0,  since  electrons  are 

bom  (ionized)  with  v * 0.  Using  Eq.  (3)  and  Maxwell's  equations  for 

h 

1 curl  E and  curl  B,  assuming  time  dependence  exp(-iuut),  and  performing 


(c2D2  + (ju2  - (Up2) (-iuo  + UD)  E = 0,  U > 0,  (4a) 

(-iiu  + UD  )(u  2(c2D2  + cu2-(U2)e=0.U<0,  (4b) 

P p ' 7 ' 


where  D = d/dz,  iup2  = 4tti»(z  )e2/my,  and  y"2  = 1-g2.  Note  that  tu  is 

related  to  the  incident  frequency  in  the  laboratory  frame  by 

(u  = (1  + 3)1^2/(1  - P)1/2,  and  u)p  is  numerically  the  same  as  the 

plasma  frequency  in  the  laboratory  frame. 

In  the  gas  (assumed  to  have  the  dielectric  properties  of  the 

vacuum),  i.e.  at  z *♦  - ",  the  solutions  are  of  course  the  incident  and 

reflected  EM  waves,  k = + tu/c.  In  the  plasma,  i.e.  at  z -*  + ®,  where 

u)  is  constant,  there  are  three  solutions:  (if  ou  2 > (U2)  the  evanescent 
P ’ P 

and  purely  growing  waves,  k = + i(u)p2  - tu2)1^2/c,  and  a propagating 
wave  k = (u/U  which  convects  with  the  plasma.  The  first  two  waves  have 
exactly  the  same  dispersion  relation  as  in  a stationary  plasma.  The 
third  wave  would  reduce  to  a time- independent  magnetic  field 
B = BQyexp(ikz),  in  a stationary  plasma,  but  in  a moving  plasma  it  is 
a genuine  transverse  EM  wave,  transporting  energy  at  the  plasma  velocity 
U.  We  shall  refer  to  it  as  the  "magnetic  wave". 

We  now  calculate  the  reflection  coefficient,  by  solving  Eqs.  (4) 
exactly.  We  recall  that  no  assumptions  have  been  made  as  to  whether 
the  front  is  sharp  or  broad.  We  consider  first  the  case  U < 0.  For 
z -*  + ®,  where  tu^  becomes  constant,  one  of  the  three  independent 
solutions  of  the  third  order  differential  equation  (4b)  is  proportional 
to  exp(ituz/U),  i.e.  represents  the  magnetic  wave.  This  solution  must 
be  excluded,  since,  for  U < 0,  it  transports  energy  in  from  z ■ + ®. 
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The  other  two  solutions  satisfy 


[d2/dz2  + (tu2  - flu  2)/c2]  E = 0,  (5) 

j 

which  is  identical  to  the  wave  equation  for  the  well-known  problem  of 
a stationary  ionization  front  in  a stationary  plasma  (i.e.  a mirror). 
The  solution  to  (5)  in  the  gas  at  z *♦  - 00  is  thus  of  the  form 

E(z)  = Ej[exp(iuiz/c)  + T exp(-ituz/c)]  . (6) 

In  general,  the  "stationary  reflection  coefficient"  T depends  on  the 
form  of  n(z),  and  the  calculation  of  T is,  in  general,  a difficult  but 
well-known  mathematical  problem.  However,  there  are  two  cases  in  which 
the  evaluation  of  |r|becomes  simple.  First,  in  the  limit  of  a sharp 
ionization  front  (i.e.  narrower  than  a wavelength), 

l-(l-<u  2/ou2)1/2 

r- E . (7) 

1 + (l-u^/flu2)1/2 

Second,  j T j = 1 if  the  plasma  is  overdense  (au^  > ou),  whatever  the 
detailed  form  of  n(z).  In  these  two  cases,  we  shall  produce  an  exact 

sn 

solution  to  our  problem.  In  the  general  case,  we  reduce  the  problem 
of  interest  to  the  calculation  of  T,  i.e.  to  the  solution  of  the 

► 

reflection  problem  at  a stationary  front  in  a stationary  plasma.  The 
[•  most  interesting  case,  of  course,  is  the  overdense  case,  in  which 

• i 

reflection  is  maximized.  We  note  that  the  criterion  for  "overdense" 


The  ratio  of  reflected  to  incident  power  (still  in  the  frame  of 
the  recombination  front)  is 


Pr/Pi  = 1 Erl 2/| Ei| 2 = |r|2.  (8) 

Lorentz  transforming  back  to  the  laboratory  frame,  where  the  plasma  is 
stationary  but  the  recombination  front  is  receding,  we  have  for  re- 
flected power 

P*/P*  - |r|2(l  - | 9l  )2/d  + I el  )2,  (9a) 

and  for  total  energy  (integrated  over  the  pulse) 


e*/e*  = |r|2(i  - |p| )/(i  + |e| ). 


(9b) 


Equations  (9a, b)  are  compatible  because  the  reflected  pulse  is  elongated 
according  to  Eq.  (2).  Equations  (5-9)  are  all  identical  to  those  for 
a receding  mirror. 

We  now  consider  the  case  of  an  oncoming  ionization  front,  U > 0. 

E is  a solution  of  Eq.  (4a)  if 


E 


exp[i(u(z-z‘')/u]g(z'), 


(10) 
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ar*;*. 


and  g is  a solution  of  Eq.  (5).  Using  the  solution  (6)  in  (10),  and 
performing  the  integration,  we  find4  in  the  gas,  i.e.  for  z - ®, 


E = Ei[exp(iojz/c ) + F [(l-0)/(l  + 0 )]exp(-iiuz/c)]  . 


Thus  the  reflected  power  is  given  by 


pr/pi  = I r | 2 ( i - 0)2/(i  + 0)2.  ( 12) 

Transforming  back  to  the  laboratory  frame,  we  find  the  power  reflection 
coefficient  is  identical  to  what  it  would  be  for  a stationary  front, 


p£/p*  * |r|2, 


( 13a) 


but  the  total  reflected  energy  is  less  than  the  incident  energy, 


e*/e*  = ! r| 2( 1 . 0 )/ ( 1 + 0). 


( 13b) 


Once  again,  Eqs.  (13a,b)  are  compatible  because  the  reflected  pulse  is 
shortened  according  to  (2).  Equations  (13)  differ  from  those  for  a 
moving  "mirror",  which  would  be 


- |r|2(i  + 0 )2/(i  - 0 )2, 


(14a) 


€r/ei 


|r|2(i  + 0 )/( 1 - 0). 


(14b) 


Our  derivation  of  Eqs.  (7-9 )>  (11-13)  has  been  quite  general; 
however,  it  is  of  pedegogical  interest  to  note  that  in  the  special  case 
of  a sharp  ionization  or  recombination  front,5  these  equations  could 


also  be  derived  in  a more  standard  way,  by  matching  gas  and  plasma 
solutions  at  the  front.  In  doing  this,  one  must  realize  that  for  the 
sharp  ionization  front  case,  there  is  an  extra  (as  compared  to  the 
usual  problem  of  reflection  at  a stationary  ionization  front)  solution 
in  the  plasma  - the  magnetic  wave.  Thus  one  more  continuity  condition 
is  needed  at  the  front,  in  addition  to  the  usual  conditions  of  continuity 
of  E and  H.  This  condition  is  that  J = 0 at  the  front  (since  electrons 
are  born  with  no  oscillating  velocity).  From  this,  it  follows  that 
&H/5z  is  also  continuous. 


We  note  from  Eq.  ( 13b)  that  the  reflected  energy  is  less  than  the 
incident  energy,  even  when  the  usual  electromagnetic  wave  in  the  plasma 
is  cut  off  (u)  < <JUp ) and  cannot  carry  away  the  excess  energy.  The 
explanation  for  this  is  as  follows.  For  U > 0,  and  in  the  limit  of  a 
sharp  ionization  front,  explicit  solution  of  Eq.  (4a)  shows5that  both 
the  evanescent  wave  and  the  magnetic  wave  are  excited  in  the  plasma. 

The  latter  transmits  energy  (both  electromagnetic  energy  and  kinetic 
energy  of  transverse  electron  currents)  into  the  plasma.  In  the 
laboratory  frame,  it  would  be  more  correct  to  say  that  magnetic 
energy  and  transverse  electron  kinetic  energy  are  left  behind  in  the 
plasma  as  the  front  moves  on.  In  the  opposite  limit  of  a broad  front. 
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kinetic  energy  of  transverse  electron  flow  convects  into  the  plasma, 
but  the  flow  of  different  electrons  is  not  coherent;  thus  the  trans- 
verse current,  and  the  amplitude  of  the  magnetic  wave,  become 
vanishingly  small.  In  either  case,  wave  energy  penetrates  an  overdense 
plasma,  and  overall  conservation  of  energy  and  momentum  can  be  demon- 
strated explicitly.5  The  external  energy  source  required  to  ionize 
the  gas  plays  no  role  in  the  wave  energy  balance. 

For  the  recombination  front  case,  U < 0,  kinetic  energy  of  trans- 
verse electron  oscillation  is  carried  toward  the  neutral  gas,  and  is 
ultimately  released  as  recombination  radiation.  Again,  it  can  be  shown 
that  energy  and  momentum  are  conserved.5  The  magnetic  wave  is  never 
excited,  since  it  would  carry  energy  in  from  z = + “. 

It  is  obvious  physically  that  EM  waves  cannot  be  reflected  from 
an  oncoming  ionization  front  moving  faster  than  c,  since  the  front 
would  instantaneously  overtake  the  wave.  What  happens  mathematically 
is  that,  in  the  super luminous  case,  the  evanescent  and  growing  waves 
in  the  plasma  both  become  (in  the  laboratory  frame)  propagating  waves. 
Both  of  these  are  excited,  rather  than  one  reflected  wave  and  one 
evanescent  wave. 

Conclusions 

Microwaves  reflected  from  an  approaching  ionization  front  will  be 
Doppler  up-shifted,  as  given  in  Eq.  (l).  The  reflection  coefficient, 
given  by  Eqs.  (13),  becomes  small  in  the  relativistic  limit  |3  1. 

Nevertheless,  it  may  be  possible  to  produce  very  impressive  power 
levels  in  the  submillimeter  regime,  using  recently  developed  high- 
power  microwave  sources.  Furthermore,  this  technique  promises  to  be 
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precisely  tuneable,  and  to  require  a very  modest  investment  in  energy 
and  equipment,  compared  to  other  methods  of  upshifting  microwaves, 
e.g.  Ref.  3. 

Finally,  we  note  that  several  calculations  appear  in  the  literature 
of  reflection  of  electromagnetic  waves  from  a moving  discontinuity 
between  stationary  dispersionless  dielectrics.6  This  model  leads  to  a 
reflection  coefficient  equal  to  that  at  a physically  moving  dielectric 
surface,  i.e.  our  Eqs.  (14),  and  thus  different  from  our  results, 

Eqs.  (12).  The  macroscopic  model  of  a dispersionless  dielectric  inter- 
face must  be  examined  carefully  to  determine  whether  it  applies  to  any 
given  physical  situation;  in  particular,  it  is  clearly  inapplicable  to 
the  case  of  an  ionization  front.  This  matter  will  be  discussed  more 
fully  in  a future  publication.5 
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Fig.  1 — The  scattering  problem  in  the  frame  of 
reference  of  the  ionization  front 
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